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Abstract 

m : 

Rectangular cavities are solvable models that nevertheless touch on many of the 
controversial or mysterious aspects of the vacuum energy of quantum fields. This 
■^j- . paper is a thorough study of the two-dimensional scalar field in a rectangle by the 

method of images, or closed classical (or optical) paths, which is exact in this case. 
For each point r and each specularly reflecting path beginning and ending at r, 
we provide formulas for all components of the stress tensor T^ u (r), for all values of 
the curvature coupling constant £ and all values of an ultraviolet cutoff parameter. 
Arbitrary combinations of Dirichlet and Neumann conditions on the four sides can 
be treated. The total energy is also investigated, path by path. These results are used 
in an attempt to clarify the physical reality of the repulsive (outward) force on the 
sides of the box predicted by calculations that neglect both boundary divergences 
and the exterior of the box. Previous authors have studied "piston" geometries 
that avoid these problems and have found the force to be attractive. We consider a 
"pistol" geometry that comes closer to the original problem of a box with a movable 
lid. We find again an attractive force, although its origin and detailed behavior are 
somewhat different from the piston case. However, the pistol (and the piston) model 
can be criticized for extending idealized boundary conditions into short distances 
where they are physically implausible. Therefore, it is of interest to see whether 
leaving the ultraviolet cutoff finite yields results that are more plausible. We then 
find that the force depends strongly on a geometrical parameter; it can be made 
repulsive, but only by forcing that parameter into the regime where the model is 
least convincing physically. 
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1 Introduction 

1.1 History and motivation 

Rectangular cavities are perhaps the most frequently studied geometries in 
connection with vacuum (Casimir) energy. Nevertheless, there are still worth- 
while things to say about them. Although exactly solvable, they exhibit many 
features that are subjects of current research and debate in the broader context 
of quantum vacuum energy: boundary divergences, corner effects, sometimes 
surprising signs, and sometimes revealing connections with geometry through 
the spectrum of periodic and other closed classical paths (or optical rays). 

A complete listing of previous literature is impossible, but we summarize what 
we see as the most important historical developments. 

Lukosz [1] calculated the interior vacuum energy of the electromagnetic field in 
a (3D) perfectly conducting parallelepiped, using zeta-function regularization, 
and predicted a repulsive (outward) force for many aspect ratios, including the 
cube. Ambj0rn and Wolfram [2] extended such calculations to a wide variety of 
dimensions, fields, and boundary conditions. Actor [5] emphasized that diver- 
gences in the total energy must be understood in terms of the local behavior of 
the field near boundaries and boundary singularities (edges and corners), and 
calculated the local zeta function for the 3D scalar field. For earlier, closely 
related discussions of a rectangular waveguide and various other systems, see 
[HIS] . All these works were done in the framework of zeta functions PTTfS] . but 
in practice, in the special case of rectangular cavities, functional equations for 
Epstein zeta functions are used to convert zeta-regularized sums over eigen- 
values into what are, in effect, zeta-regularized sums over classical paths. (In 
[H] a transition to an ultraviolet cutoff was also made at this step.) 

The calculation of vacuum energy via classical paths (also called optical rays) 
|lU|lllfl2[13|ll4fl5lll6fl7] leads naturally to more physical regularizations as- 
sociated with separation of points in the Green functions of the theory. For 
rectangular parallelepipeds such calculations are exact (no stationary-phase 
approximations are required) and reduce to the classic method of images. Al- 
though they did not discuss vacuum energy, Balian and Bloch [TS] used the 
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3D parallelepiped as a principal example in their landmark study of the re- 
lation between periodic orbits and oscillations in the eigenvalue density, and 
their catalog of periodic and other closed orbits is the best starting point for 
a study of the rectangle. Hacyan et al. [12] calculated the full stress tensor 
for the electromagnetic field in a box by a Green-function approach; along 
the way, they showed how to reduce the electromagnetic problem to scalar 
fields with mixed Dirichlet and Neumann boundary conditions via Hertz po- 
tentials. In this connection see also [20]. (Contrary to the impression left by 
some papers on the global problem, a local investigation of electromagnetism 
cannot be split into pure Dirichlet and pure Neumann problems, even when a 
decomposition into TE and TM modes exists.) 

At this point we should mention the work of Ford and Svaiter [21], which 
showed that physically motivated cutoffs could convert divergences into finite 
effects clearly localized near boundaries. This theme has been repeatedly vis- 
ited since then [2"2"f2"3][T5 ,24,25,26] and will play a major role in the present 
work. 

Cavalcanti [27] rejuvenated the field by introducing the piston model (for a 
2D scalar field), discussed in detail in Sec. Eland illustrated in Fig El (Similar 
ideas were advanced earlier by Svaiter et al. [28].) The motivation for the 
piston is that the calculation of the force on the piston plate is unaffected by 
either uncompensated divergences or unknown forces from the exterior. The 
conclusion of [2Z] is that the force is always attractive (inward). That paper 
used both zeta and cutoff regularization, but still starting from the eigenvalue 
spectrum. Hertzberg et al. [291130] extended the piston model to dimension 3 
and to the electromagnetic field, and they analyzed it in terms of closed paths 
(but without the close attention to locally defined quantities that we provide 
here). From that point of view, the repulsive nature of the Lukosz force is 
attributable to a particular type of path moving parallel to the plate and 
producing an energy proportional to the piston displacement. (It is essentially 
the Casimir energy associated with the walls perpendicular to the movable 
plate.) But such energy is also present in the exterior part of the piston shaft, 
and therefore these paths make no net contribution to the force. What is left 
of the Lukosz force is attractive. This effect shows up even more clearly in the 
two-dimensional model (Sec. H]). 

Rodriguez et al. [3~Tf32] have made a numerical study of two conducting rect- 
angular objects in a narrow passage, a model closely akin to the pistons and 
pistols we discuss here. They conclude that the distance to the confining walls 
influences the attraction between the blocks, and their analysis makes use of 
the local stress tensor. In [33J that model is approached by the method of 
closed optical paths. 

Illuminating though the piston has been, it does not settle the original issue 
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of the physical reality of the force calculated by Lukosz pQ and others. The 
existence of a Casimir-like energy in the exterior part of the piston shaft says 
nothing about what happens when that part of the shaft is removed, the plate 
remaining free to move (see Fig. [7|in Sec. The "finite part" of that force 
is robust, in the sense that all reasonable prescriptions for calculating it give 
the same answer. It can be obtained by differentiating the total energy, or 
by integrating the pressure over the movable boundary. It can be obtained 
by zeta-function regularization or by ultraviolet cutoffs, and within the latter 
framework the choice of cutoff function dictates the relative sizes of the cutoff- 
dependent terms but not the structure of the series nor the numerical value of 
the finite term [27| Appendix B]. Is the object of this consensus a meaningless 
number? One of our goals is to investigate to what extent it has physical 
significance. 

The opinion expressed in [22] is that "Without [the piston shaft] (or some 
open region that allows rigid motion of the partition) the Casimir energy of the 
parallelepiped is, in fact, cutoff dependent. If the cutoff dependence is somehow 
ignored, a repulsive force . . . remains as an artifact." We agree that a correct 
calculation of the force on the piston must include the effect of the piston 
shaft, and that the net effect is attractive. We do not agree that the repulsive 
force associated with the interior can be dismissed as an artifact of naively 
discarding a divergent term. The scenario indicated in Fig. [7J a box with a 
movable lid, is a well-defined problem of relative motion of rigid bodies, just 
as much as the piston is; cutoff-dependent energies associated with the rigid 
boundaries cannot affect the force. The difficulties of analyzing Fig. Hare, first, 
that the effects of the corners and gaps in that geometry are hard to calculate 
(see, however, [34 35 36 |37|l38j and, on a different tack, [391I4TJ] ). and, second, 
that the idealized Casimir theory is not physically applicable to very small 
separations of the bodies. We evade the first problem by considering another 
scenario, the "pistol" (Fig. [Sj), which should still exhibit the uncompensated 
Lukosz force on a flat boundary. However, we find that the situation is then 
confounded by a strong countervailing attractive force associated with the 
Casimir energy in the narrow gap surrounding the "bullet". Moreover, one now 
runs up against the second problem, which cannot be treated seriously within 
the limits of our methodology. A crude model of a "real" boundary can be 
easily obtained, however, by maintaining a finite cutoff of atomic dimensions. 
The result is that the force depends sensitively on how tightly the "bullet" 
fits into the "barrel". If the fit is loose, the Lukosz force is overwhelmed by the 
corresponding force associated with the gap surrounding the bullet, and the 
net force is attractive. If the fit is tight, the gap force can be made repulsive 
or even fine-tuned to vanish, as originally hoped; unfortunately, that is the 
regime in which one is least justified in taking the model seriously. All we 
claim is that external forces opposing the Lukosz force are model-dependent 
and might, in principle, be controlled so as to demonstrate the existence of 
the Lukosz force. 
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In this paper we consider strictly the two-dimensional scalar field, usually 
with Dirichlet boundary conditions, although we sometimes lapse into the 
three-dimensional electromagnetic terminology (such as "conductor" ) for con- 
ceptual discussions. It is intended that three-dimensional generalizations will 
be presented elsewhere [H] . Sec. [2] presents a thorough analysis, by means of 
classical paths, of all components of the stress (energy-momentum) tensor in 
a rectangle. Sec. [3] does the same for the energy and also the pressure and 
force on one side. Contributions are recorded for each path (or class of similar 
paths) separately, with comments on their physical or geometrical significance 
to the extent that we can discern it. The results are stated for all values of the 
curvature, or conformal, coupling constant, £ (see (j5J) and and all values 
of the parameter in an exponential ultraviolet cutoff. For the most part, they 
are stated for any combination of Dirichlet and Neumann conditions on the 
four sides of the box. A brief account of this part of the work was published 
previously [12], along with evidence that the gravitational effects of bound- 
aries in the "renormalized" theory without cutoff can be understood (and 
believed) as the distributional limit of the predictions of the cutoff theory, 
thereby providing a true renormalization. In the rest of the paper we restrict 
to the Dirichlet condition. The piston is reviewed from our point of view in 
Sec. |U Sec. introduces the pistol model and treats it in the Dirichlet theory. 
Sec. [6] investigates the pistol with a finite cutoff. 



A remark on terminology: Many authors, including some of ourselves on pre- 
vious occasions (e.g., [23]), use the term "renormalized energy" to refer to the 
finite part of a regularized energy when the latter is expanded as a series in the 
cutoff parameter. Strictly speaking, "renormalization" refers to the process of 
obtaining physically observable probability amplitudes by absorbing suitable 
divergent and finite contributions into redefinitions of physical parameters 
(couplings, masses, etc.) appearing in the bare Lagrangian. Ideally, all renor- 
malizations in the first sense should either be associated with renormalizations 
in the second sense or be justified by cancellations of divergent terms coming 
from different sources. Yet in the absence of a completed theory, one must 
often talk about renormalization in the first sense without having an obvious 
counterterm or cancellation, and there seems to be no convenient substitute 
terminology. Much of our work in this paper has to do, in fact, with exhibiting 
cancellations, and [12] and its planned sequels have to do with gravitational 
counterterms. When we use "renormalization" or "renormalized" in the first 
sense, we have always either put the word in quotation marks or accompanied 
it by the word "naive", depending on context. 
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1.2 Basic formalism 



We are concerned here with the massless scalar wave equation 

d{x°y v 



in a cavity Q together with a Dirichlet (<fi = 0) or Neumann (n ■ V0 = 0) con- 
dition on each part of the boundary of Q. We write H for the corresponding 
positive self-adjoint operator: H = — V 2 with boundary conditions under- 
stood. The eigenvalues of H are positive, with the possible exception (in the 
totally Neumann case) of a constant eigenfunction with eigenvalue zero. The 
formulas in this subsection are presented for arbitrary spatial dimension d, 
but in the next section we specialize to d — 2. 



The field equation 
and Lagrangian 



JT]) is obtained canonically from the curved-space action 



S 



L 



d d+l : 



(2) 



by taking the variation with respect to <p and then setting the metric to its 
flat-space value, g^ u = rj^ . (Our tensorial sign conventions are that r] 0Q < 0, 
but T o > for normal matter.) The stress tensor is defined by 

r = ~ (3) 
" 1 og^v 



It reduces in flat space-time (after use of the equation of motion, ([I])) to 

T^u = (1 - 20 d v <f> + (2f - \) Via , d x <j) d x <J> -2^<j> d^d u <p. (4) 



In (J2J) R is the curvature scalar, and £ labels different possible gravitational 
couplings. In curved space different values of £ are different theories; after the 
reduction to flat space the field equation is independent of £, but the stress 
tensors are different. It turns out [see ([7j)] that changing £ changes T o only 
by a divergence, and therefore the total energy E = J Q T 00 dr is independent 
of £, at least classically, under Dirichlet or Neumann boundary conditions. (A 
Robin boundary condition [43. 23J44] , n • V0 = 70, would require a boundary 
term to be added to the action There results a ^-dependent boundary 
term in E, which vanishes when £ = | . Similar remarks apply to models with 
delta function potentials j24~f25f26j.) There are three natural choices of £: 
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£ = : minimal coupling, which simplifies the Lagrangian and curved-space 
field equation; 

£ — £d : conformal coupling, 



d-l 
Ad 



^2 a j £3 k ) £00 4 j 



(5) 



which results in the mildest behavior of the quantized field near the bound- 
ary; 

£ = \ : the coupling that eliminates the Robin boundary energy, which also 
simplifies the relation between the stress tensor and the total energy, as we 
shall see. 

It is convenient to adopt £ = \ as the base value and to define (5 — £ — \ to 
parametrize the coupling. Thus we write 

T M ,(0 = ^(i) + AT^ (6) 
and obtain 



T 00 ( 



\dx° J 
d(p d<f> 



dx° dxj 



0V 2 



d 2 



AT oo = -2/3V-(0V0), 



dx° dx-i 



0.7 



-2(3 



8 ( d<f> 
dx° \ dxj 



(7) 



, (8) 



dcj) d(p 



dxj dxk 



AT jk = -2p 



d 2 <j> 



dxj dxk 



+ </>■ 



d 2 



dxj dxk dxj dxk 



when j k, 



(9) 



'dcj) 
dxj 



AT j3 = -2? 



d<j> 



dx° J 

The trace of the tensor is 

n = - (i + w) 



<9 2 0" 

dxj 2 



E 



' d<j> 
dx k 



+ 



d 2 f 
dxj 2 



(10) 



dx° 



- (V0) ; 
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which vanishes for the conformal coupling, (3 = — (4d) 1 . 

When the theory is canonically quantized, the vacuum expectation value of 
the stress tensor is expressed formally in terms of the normal modes 

ip n = -jj= n (r)e- iuJ " x ° , Hcj) n = u£(j> n , \\(j) n \\ = 1, 

V 



as 



&(r)> = E^^n,<]- 



n=l 



'121 



(The Neumann zero mode, when it exists, is omitted and ignored. If included 
and treated properly, it would add a continuous energy spectrum [IS, Ap- 
pendix].) The notation in (Tl2l) means that in each of the bilinear terms in 
fJ7|)- f|T0|) . the field is replaced by a mode function in one factor and by 
its complex conjugate in the other. (When the factors are not the same, the 
product should be symmetrized.) In particular, 

(Too(|)) = iE^|0n(r)| 2 . (13) 



Integrating T o(|) over Q gives the expected formal sum for the total energy, 



As promised earlier, we regularize all these divergent sums with an exponential 
ultraviolet cutoff. It is convenient to start from the (Poisson) cylinder kernel, 

T(t, r, r') = £ n (r)0 n (r')*e-^ = (r| e -^|r'>. (15) 

n=l 



(Here t is not the physical time.) Then 
1 dT 

(r o(|)>t = -2^(i,r,r), (16) 
(E) t = -\ ^T(t), T(t) ee / T(t, r, r) dr. (17) 
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To obtain (AToo) and the other components of (T^) one needs a more prim- 
itive cylinder kernel, 



T(t, r, r') = - £ 



n=l U n 



Then T = ^ and 



(AT 00 ) t = i3V F -[V rl T{t,r,r') 



(19) 



In terms of partial differential equations, T and T are characterized by the 
elliptic equation 



d 2 T 



-V 2 T 



(20) 



along with the imposed spatial boundary conditions, the initial condition 
T(0,r,r') = 5(r-r') = ^(0,r,r'), 



and the requirement of boundedness as t — > +00. (The Green function T 
can be introduced differently, either as twice the Euclidean Green function 
in R x Q with its source on t = 0, or through an analytic continuation to 
imaginary time of the Wightman or Feynman two-point function.) 

The vacuum expectation values of the summands in (jSJ) are identically zero, 
as expected from the mode-by-mode time-reversal invariance. For the other 
components one obtains 



(T 3] {\))t 



-2 



d 2 



d 2 d 2 

+ TT^ + 



2 



dxj dx'j dxj 2 dx' 2 
2 



T. 



d 2 



d 2 d 2 

+ 7^ o + 



dt 2 dxk dx' k dxj 2 dx'f 

d 2 d 2 d 2 



+ 



T. 



d 2 



8 [dxj dxk dx'j dx' k dxj dx' k dx'j dx k 



T. 



Q2 



+ 



d 2 



+ 



d 2 



+ 



d 2 



dxj dxk dx'j dx' k dxj dx' k dx\ dx^ 



T. 



(21) 
(22) 

(23) 
(24) 



where it is understood that r' is to be set equal to r at the final step. 
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2 The stress tensor 



2. 1 Preliminaries 



We now restrict attention to dimension 2 and write x for x\ and y for x 2 • 
Define 



^f = S, (25) 



^ _dt 2 

_ 1 /<9 2 T <9 2 T\ _ 1 (d 2 T d 2 T\ 

d 2 T „ ,9 2 T 



^ dxdx 1 ' ^ 2 dydy'' 

_ i / a 2 T a 2 T \ 

12 2 \ <9x <9y' <9y dx' J ' 

__ 1 / d 2 T d 2 T \ 
12 2 \dx dy dx' dy' J 

(The subscripts on D and E are merely to facilitate later generalization to 
higher dimensions.) Then from f|T6|) . (fT9l) . and ([2T]) - (|24p we have 



(T 00 (i)) t = -|A (30) 

(AT 00 ) t = /3(S! + 5 2 + d + C 2 ), (31) 

(r ol (l)> 4 = 0=(Ar ol > tJ etc., (32) 

(TnmMiBi-Cil etc., (33) 

(AT n ) t = (3(A + B 1 -C 2 ), etc., (34) 

(T 12 (i)), = i(-A 2 + £i 2 ), (35) 

(AT 12 ) t = /3(D 12 + E 12 ). (36) 

In ( I30|) -(l36l) it is understood that r' = r. 

2.2 Path classes and energy density 

The cylinder kernels in infinite two-dimensional space are 

T(t,r,r') = -^(t 2 + |r-r'| 2 )- 1 / 2 , (37) 

Z7T 

T(t,r,r') = f (t 2 + |r-rf)- 3 / 2 . (38) 
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x = point r under study, 

• = periodically displaced image, 
o = reflection through a side, 

* = reflection through a corner. 



Fig. 1. A point in a rectangle and its images relevant to Dirichlet and Neumann 
boundary conditions (cf. |181 Sec. 9.A]). Image points fall into three classes according 
to whether the number of reflections is even in both dimensions, one, or neither. 
The first case corresponds to periodic displacements. Points of the third class are 
joined to r by lines that pass through an intersection point of the lattice of extended 
rectangle sides — i.e., an image of a corner of the rectangle. 





T 










i 
i 














• 















Fig. 2. Two periodic paths (one solid, one dashed) are shown, both within the 
rectangle and in the covering space. 

Because of its central importance, we shall discuss the energy density, (T 00 (|)}, 
along with the construction of the cylinder kernel for the rectangle, path by 
path. For rectangular parallelepipeds of any dimension, with any combination 
of Dirichlet, Neumann, and periodic boundary conditions, the construction of 
any kernel (Green function) as a sum over classical paths reduces to the classic 
"method of images" (and yields the exact answer). For a rectangle the array 
of image points appears in Fig. [TJ To every path is associated a sign, (—l) v , 
where r] is the number of Dirichlet sides struck by the path. (If a path hits a 
corner, both sides are counted, and the path reflects back upon itself.) The 
image sum for T is not absolutely convergent, but the derivatives of the series, 
from which observable quantities are calculated, are convergent. 

Following Cavalcanti [27] we take the rectangle to have horizontal and vertical 
dimensions a and b, horizontal and vertical coordinates x and y, and horizontal 
image-displacement indices j and k. (We occasionally still find it necessary to 
use j and k as tensor indices, but never in the same equation as the image 
indices.) Thus the contribution of a typical periodic path (see Fig. [2]) to T is 

Tp jk = ~ + (2ja + x'- xf + (2kb + y'- y) 2 )}- 1/2 . (39) 



x 



From (1301 and (1251) we obtain 
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Fig. 3. Two corner paths. Inside the rectangle such paths bounce back from a corner 
and retrace themselves. The shortest such paths have lengths arbitrarily close to 0. 

(-1)" 



tPjk 



An 



t 2 + (2jay + (2kbY 



-5/2 



-2t 2 + (2ja) 2 + (2kb) 2 
which is independent of r. Also, one finds from (13T1) and 



(40) 

-(1271) that Cj = 



-Bj in this case and hence 



(AT 00 ) 



tPjk 



0. 



(41) 



These two results are expected and related: Since AT 00 is a total divergence 
and hence must integrate to 0, and since the energy from a periodic path is 
independent of position in the rectangle, (AT o) tP - fe must be identically zero. 
Finally, note that if the boundaries are all Dirichlet or all Neumann, rj is even 
and hence (Too)opjfc * s always negative. 



Next, consider the simple and interesting case of a corner path: 



2tt 



t 2 + (2ja - x' - x) 2 + {2kb - y' - y f 



-1/2 



(42) 



This time one finds that C,- 



-Bj and B\ + B 2 = —A, so that 



CT«,(i)> 



tCjk ~ 4/3 (^Too) t Cjk 

1 k 2 + (2ja - 2x) 2 + (2kb - 2y) 2 ] ^ 



x 



47T L 

~-2t 2 + (2ja - 2x) 2 + (2kb - 2y f 



(43) 



That is, the two terms in (7oo) tC - fe are proportional, and, in particular, 
(Too) tC j k vanishes for minimal coupling (/3 = —\)- These seeming coincidences 
are probably related to the fact that the integral of (T 00 ) tC j k over the rectangle 
must vanish (see Sec. |3]). Note that the quantity is a function of the distance to 
r from the corner or corner- image concerned (see Fig. [3]). Again it is negative 
as t — > whenever all the boundary conditions are of the same type. 



The situation is slightly more complicated for paths that "bounce" in one di- 
mension while being periodic (or fixed) in the other. The number of reflections 



12 







o 




o 




? 

1 












o 




o 




o 




o 





Fig. 4. Two side paths of the vertical subclass. The dashed path is a direct perpen- 
dicular reflection (retracing itself), with length approaching as r approaches the 
top boundary. The solid path combines a reflection from the top with a periodic 
horizontal drift; its length is bounded away from 0. 

is now odd, and the energy density turns out to be positive for Dirichlet con- 
ditions and negative for Neumann, at least for (3 near 0. We call these paths 
"vertical side paths" if the bounce is off a horizontal side (see Fig. HJ); this 
includes, in particular, the strictly vertical paths (j = 0). In those cases we 
have 



T 



Vjk 



2vr 



t 2 + (2ja + x' - xf + {2kb - y' - yf 



-1/2 



(44) 



■l) v r 



tVjk 



47T 

~-2r 



t 2 + (2ja) 2 + (2kb - 2yf 
f (2ja) 2 + (2kb - 2y) 2 



-5/2 



(45) 



B(— IV? r 

(AToo W = [t 2 + (2ja) 2 + (2kb - 2y) 

t 2 + (2ja) 2 - 2{2kb - 2yf 



-5/2 



(46) 



These quantities depend only on y, not x\ in other words, such a term is a 
function of the distance from a wall or an image of a wall. In this case the two 
terms in the energy density are distinctly different, so it pays to write out the 
total explicitly: 



(Too) 



:-D" r 



tVjk ' 



7T 



X 



t 2 + (2ja) 2 + (2kb - 2yf 
(2ja) 2 - 



-5/2 



(2/3 + J) (2kb - 2y) 2 ] . (47) 



The most interesting observation here is that the coefficient of [2kb — 2y) 2 
vanishes for conformal coupling (/3 = — |). When j = and k = or 1, the 
energy density for t = generically has O (y~ 3 ) divergences at the boundary, 
but those divergences are removed in the conformal case; this is as close as one 
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comes in a rectangle to the well known fact that the energy density between 
infinite parallel plates is constant in the case of conformal coupling. 

Formulas for horizontal side paths are easily obtained from (I44 p - (|47l) by in- 
terchanging the roles of the two dimensions. 



2.3 The other components 



From f[3"51) - f[3"oT) and f[3"9"j) . fT4"2"j) . (jHj), one finds the spatial components (pressure 
and shear stress) . We omit the formula for the 22 component when it is obvious 
from the 11 formula. 



Periodic paths: 



(TiiU))tPjk 



(Ti 2 {\))tPjk 



-IV r 



47T 



t 2 + (2ja) 2 + (2kbf 



-5/2 



t 2 - 2(2 3 a) 2 + (2kb) 2 
3(-lY> 



71 



t 2 + (2jay + (2kb) 2 jakb, 



-5/2 



(ATn)tPjjfc — — (ATi 2 )tpjk ■ 



(48) 

(49) 
(50) 



Thus the stress tensor associated with a periodic path does not depend upon 
the conformal parameter, nor upon the coordinates. The individual terms (T12) 
are nonzero, but they add to zero when summed over either j or k, as reflection 
symmetry requires. 

Corner paths: 

(TiiCi))*** = (T 12 (i)) tCife = 0, (51) 



<AT n > 



tCjk 



71 



t 2 + (2ja - 2xf + {2kb - 2y f 



-5/2 



t 2 + (2ja - 2x) 2 - 2(2kb - 2yf 



(52) 



(AT 12 ) 



tCjk ■ 



7T 



x (ja — x)(A;6 — y). 



(53) 
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In addition to (and in contrast to) the remarks about the energy density 
made below (l4~3]) . we observe: (1) The spatial components of the corner-path 
stress tensor vanish when £ = | (whereas the energy density vanishes when 
£ = 0). So far we have no intuitive explanation of this fact. (2) The spatial 
components are no longer functions of corner-image distances alone, though 
they do have (for t = 0) an 0(|r|~ 3 ) dependence on corner- image coordinates, 
as the energy density does. (3) When (3 ^ there is a nonzero (T12), which 
does not vanish even when summed. However, if we evaluate it on a boundary 
(such as x = (integer) x a), where it would have a clear physical interpretation 
as a shear force on the wall of the box, then it does vanish when summed. 



Vertical paths: 



(Tn(j))tvjk 



-!)*> 



t 2 + (2ja) 2 + (2kb - 2yf 
2(2ja) 2 + {2kb - 2y) 2 



-5/2 



(54) 



(AT n ) 



tVjk 



P{ l)V \t 2 + (2ja) 2 + (2kb-2yf' r ' 2 



x 



7T 

t 2 + (2ja) 2 - 2(2kb - 2y) 2 
\T2i)tvjk = = (Ti 2 )ivjfc = 0. 



(55) 
(56) 



In addition to the remarks surrounding fT4"4"|) - (HT|) . observe that (T u2 ) = 
for all v. That is understandable: there is otherwise no way to satisfy the 
conservation laws fl60l) for \i = 1 and fi = 2 by functions that depend only on 
y but are not constant. 

Horizontal paths: 

(Tn)m jk = = (T 12 ) mjk = 0, (57) 



(T 2 2{\))mjk 



4)17 



47T 



t 2 + (2ja - 2xf + (2kbf 



-5/2 



x 



t 2 + (2ja - 2xf - 2(2A;6) 2 



(5f 



(AT 22 )tHjfc — — 



7T 



t 2 + (2ja - 2x) 2 + (2A;6) ; 



-5/2 



x 



t 2 - 2(2ja - 2x) 2 + (2kbf 



(59) 
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Observe that (T12) = for all side paths. 

For the formulas above one can verify the conservation law 

- + h Ti ' ] + = ( " =0,1 ' 2) - (60) 

Here the first term is always 0, because the quantities do not depend upon 
time (not to be confused with the regularization parameter t). In the conformal 
case, P = — I , one also has tracelessness, 

- (T 00 ) + (Tn> + (T 22 ) = 0. (61) 
These identities hold for all t, not just t — 0. 



3 Energy and force 



3.1 Introductory remarks 



In this section the results of the previous one will be used to calculate the 
contribution of each image term to the total energy, E, of the scalar field in 
the rectangle, and consequently the force, — on the rectangle's right side. 
We are concerned here only with the force from inside the rectangle; "piston" 
arrangements in which it is possible to calculate or estimate forces from outside 
will be considered in later sections. Consequently, uncompensated divergent 
terms arise as the cutoff parameter t is taken to 0, and such terms need to be 
identified and systematically isolated for later physical scrutiny. 

The sign of Casimir energies and forces has long been a topic of great interest 
and mystery, and one of the motivations of our research has been to see what 
light the decomposition into image terms, for which the sign is easy to under- 
stand, can shed on such questions. The following discussion is easy to present 
for arbitrary spatial dimension d. 

The cylinder kernel in R d is 

T(t, r, r') = C(d) t(t 2 + |r - rf )-( rf +D/ 2 , C(d) = -^2 . (62) 

7T 2 



Consequently, in the ci-dimensional analogues of the constructions in Sec. 12.21 
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all the terms in the energy density (TboQ)) will have the form 



2 dt L " (t 2 + W) s+1 



where W is some nonnegative function of r. If W > 0, the limit as t — > is 
— t;(—1) v W~ s ; for pure Neumann boundary conditions it is always negative, 
while for pure Dirichlet conditions it will be positive whenever the number 
of nonperiodic (bounce) dimensions is odd for the path concerned. If W — 0, 
the small-t behavior is (—l)' n (s — |)t~ 2s , divergent and opposite in sign to the 
other terms. 

Now we consider integrating over a coordinate u when W = V + {mL — u) 2 
with V > 0: 

L r {s-\)t 2 -\[V + {mL-uf] 
J [t 2 + V + {mL- w) 2 p+! 



By the mean value theorem for integrals, 

As-\)t 2 -\[V+{mL-Cf] 
[t 2 + V + {mL - C) 2 ] s+1 ' 



where < £ < L and £ may depend on £. But for us s is always greater than 
or equal to 1. So if the integral converges at all, the integrand is bounded and 
we are back to the situation of the previous paragraph with W{t) > and 
having a positive lower bound. In the contrary case, V — and m — or 1, 
the situation is more delicate but the sign question is absorbed into the issue 
of the physical meaning of surface divergences. 

In summary, for all finite terms we have good control over the sign. (The 
(AT o) terms are irrelevant to total energy, as discussed below.) One under- 
stands why Ambj0rn and Wolfram [21 Table I] found nontrivial signs only 
for Dirichlet problems (not periodic or Neumann), and the particular sign 
patterns they saw are not surprising. 

To understand the significance of various paths, it is useful to refine the clas- 
sification of paths in the previous section. Each closed path is characterized 
by its image indices, j and k, and by its periodicity type, P, V, H, or C. 

P: Periodic paths, producing constant terms in the energy density 

PZ: j = = k — the zero-length path 
PV: j = 0, k 7^ — vertical periodic paths 
PH: k = 0, j 7^ — horizontal periodic paths 
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PD: j ^ 0, k ^ — diagonal periodic paths 

V: Nonperiodic closed paths whose uncompensated "bounce" occurs on the 
top or bottom side of the rectangle, producing energy densities depending 
on y only 

VP: j = — perpendicular vertical bounce paths 

VD: j ^ — vertical bounce paths with horizontal periodic drift 

H: Nonperiodic closed paths whose uncompensated "bounce" occurs on the 
right or left side of the rectangle, producing energy densities depending on 
x only 

HP: k = — perpendicular horizontal bounce paths 

HD: k 7^ — horizontal bounce paths with vertical periodic drift 

C: Closed paths that are periodic in neither dimension, producing energy 
densities associated with corner images 

Path PZ produces, by (ITT)) , the ubiquitous volume (here area) divergence, 

^ / x ah , ab , . 

which, being ubiquitous, is always ignored (except for possible relevance to 
cosmological dark energy). All other terms in the energy density are pointwise 
finite, but some of them have nonintegrable divergences at the boundary. The 
path classes involved are VP and HP, which produce the well known surface 
(here perimeter) divergence in the total energy, and C, which produces an 
energy density that seemingly diverges at the corners but nevertheless makes 
no contribution to the "renormalized" total energy, as we shall see. 



3.2 Energy calculations 

Let us first dispose of fj (AT 00 ) dxdy, which is expected to be zero because 
AT o is the divergence of a vector field, proportional to 0V0, that vanishes 
on every Dirichlet or Neumann boundary. From (|41j) and (TT3)) we see that 
the quantity is indeed zero for periodic paths, while for corner paths it is 
proportional to the T 00 (±) term (which also will turn out to be zero). The 
situation for side paths is more subtle. The integral of (AT o) tVj - fe from (T46T) 
is not zero, which is not surprising since the field from a single image source 
does not satisfy the boundary conditions. However, because (TT6j) is a total 
derivative, a calculation almost identical to that in (IT9j) - (l80)) below shows 
that the sum over k does telescope to 0, at least when the top and bottom 
boundaries are of the same type (both Dirichlet or both Neumann). 
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The total energy contributed by a periodic path is trivially obtained by mul- 
tiplying (flOl) by the area, ab. The sum of all such terms splits into PV, PH, 
and PD parts as 



ab » (2fc&) 2 - 2t 2 a6 ~ (2ja) 2 -2t 2 



27r fri [t 2 + (2A;6) 2 ] 5 /2 27r^ v ' [t 2 + (2ja) 2 ] 5 /2 

a6 vvr IV (2j^) 2 + (2fc&) 2 - 2t 2 

vr ^ fc tf ] [i 2 + (2ja) 2 + (2kbff/ 2 ' ^ ^ 

If all four sides are of the same type, ( 1651) simplifies in the limit t — > to 

v 7 x 16tt6 2 167ra 2 8tt f-f \ J > K ' 



(a well known result — e.g., [7]). In the special case of a square of side a, 
numerical evaluation of fl66l) gives — 0.089859/a (identifiable with the vacuum 
energy of a torus of dimension 2a as recorded in [2|7] ). Because of the need to 
sum over a two-dimensional lattice, the numerical convergence is rather slow, 
even when repetitions of primitive orbits are handled all at once — in contrast 
with the situation for parallel plates, where the sum over paths has been found 
to be very efficient and increasingly so in higher dimensions [T6P6] . 



To calculate the other Too(f) terms it is convenient to return to the cylinder 
kernel T and integrate it over the rectangle before taking the final t derivative 
in d25]). 

Consider first the paths of subclass VP. According to (H4|) with r' = r, for any 
path of class V we have 



T V j fc (t,r,r) = + 



2tt 



t 2 + (2ja) 2 + (2kb - 2yf 



-3/2 



(67) 



Setting j = we arrive at 



a b 

j dx J dyT Yok (t, r,r) 
o o 



2tt 



at 



-i — 3/2 

t 2 + {2kb-2y) 2 dy. 



The terms with k = and k — 1 are divergent (when t — > 0) at the bottom 
and top boundaries, respectively. The other cases are finite, but need to be 
added to the divergent ones to build up a "clean" divergence, proportional to a 
power of t, that can be discarded in a systematic renormalization of the mass of 
the boundary plate. (Formally, this quantity is the total energy of an isolated 
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surface in otherwise empty space [3].) Again we consider for simplicity only 
the case where both horizontal boundaries are the same type, so that (— l) v 
is independent of A;. As a reminder that this assumption is in force, we shall 
write in resulting formulas 



The terms combine easily (telescope): 



b 



E JdxJ dyT Y0k {t, r, r) = / (t 2 + Ay 2 ) 3/2 d y = T ^- . 



k=—oo 







Obviously the formula for class HP is the same with a replaced by b. Therefore, 
the total contribution from VP and HP to the trace of the cylinder kernel can 
be written as 

TxM («» 



where P is the perimeter of the rectangle. It corresponds to a divergent surface 
energy 

W" = =FgSS- (™) 

For paths of class VD we obtain 

00 r } at I t 

E y & y d V T Vjk (t, r, r ) = T - — = , + o 



fc=— oo 







2tt t 2 + (2ja) 2 87raj 2 \f J ' 

The sum over j gives the well known C(2), so 

oo a b , 
EE J dxj dyT Yjk (t,r,v) = T ^ + 0(t 3 ). 

j^0k=~oo 

The corresponding contribution to the energy is ±7r/48a; it may be thought 
of as a Casimir correction to the surface energy of the sides at y — and y = b 
caused by the presence of the perpendicular sides with separation a. Thus the 
total energy from VD and HD paths is (at t — 0) 

(E)tB = ±^(- + \). (71) 



48 \a b 
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It is comparable in magnitude to the term from periodic paths, fl66l) . In fact, 
for the square it is larger, since 7r/24 m 0.13; that is why the "renormalized" 
(Lukosz) energy of the Dirichlet square comes out positive. (The situation for 
the force is different, however, as we shall see.) 



Finally, for a path of class C we have from 



T Cj h{t, r,r) 



-1)" 
2tt 



t 2 + (2ja - 2x) 2 + (2kb - 2yf 



-3/2 



(72) 



Terms with {j, k} C {0, 1} yield divergent integrals in the energy ( ~ / r~ 3 r dr) 
if t is set equal to 0, but if one integrates with t positive, the result is quite 
different. We assume that all sides are of the same type, so that (— 1) T? = +1. 
Then the contribution to the cylinder trace from the corner paths telescopes 
to 




Being independent of t, this term makes no contribution at all to the energy 
via (l30j) . (In a related independent calculation by Zaheer et al. [33] the corner 
paths were not even considered, because the rectangle was obtained as a lim- 
iting case of a configuration where they did not exist.) In the next subsection 
we shall review why this result is exactly what should have been expected. 



3.3 Relation to heat kernel asymptotics 



Let K(t, r, r') be the heat kernel corresponding to the system under study 
({r\e~ tH \r') in quantum-mechanical notation, as contrasted with (|T5|) ). Let 
K(t) be its trace (cf. (Il?p ). It is well known [4"T|48II49|I8] that as t — > 



A P 1 

Kit) = — — =F — 7= + - + 0(t c 



(74) 



where A = ab and P = 2(a+b) are the area and perimeter of the rectangle, and 
=1= is as in (|68l) . (Here we state (fT4"]l only for the cases where all four sides are 
of the same type. The other cases — in which, for instance, the second term 
is not proportional to P, but the qualitative conclusions of this subsection 
remain true — are discussed in [50].) 
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It follows [23|[5Tj that the trace of the cylinder kernel must have the expansion 

and hence by (fTTj) the regularized Casimir energy is 
1 f)T A P 

where E ren is a constant traditionally identified as the "renormalized" Casimir 
energy. (E reri is not determined by the heat kernel expansion (1T41 : it is hidden 
in the 0(t°°) term there.) 

Our calculations above have confirmed ( 1751) and ( 1761) and determined E ren . 
The terms in ATS]) are AMI). fl55} . and (1751). £ ren is the sum of flSSD and <^> 
(in agreement with previous authors, including [2,27*33]). 

With regard to the inevitability of the disappearance of the corner energy 
(without an explicit renormalization of any kind!), we stress [23|5T] that the 
coefficient of 1/t in (1761) must be 0. (For dimensional reasons, that is where a 
corner term would need to appear, along with contributions linear in boundary 
curvature or in a Robin constant.) A f 1 term in E(t) would have to come 
from a In if: term in T(t), which in turn would be associated with a logarithmic 
term in K(t), and such terms do not exist. On the other hand, there is no 
general reason why E(t) could not contain a hit term (and a resulting scale 
ambiguity in the "renormalization"). That would correspond to a tint in T(t) 
and hence a t 1//2 in K{t) — which can actually occur (for example in a disk), 
but does not in the model under study here. 



3.4 Force and pressure calculations 



We now investigate the force on the side at x = a from the field vacuum inside 
the rectangle, in the case where all sides are Dirichlet. From the previous 
subsections, the naively renormalized energy yields the force 

ggren Cg) C(3)6 & f» k 2 b 2 - 2jV 7T 

da 16nb 2 87ra 3 8tt (fa 2 + k 2 b 2 f' 2 48a 2 ' 1 ' 



where the terms are the contributions of path classes PV, PH, PD, VD, and 
HD, respectively. 
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It is important to remember that positive energy is not always the same thing 
as positive (repulsive) force, although that is true in many of the classic 
Casimir-force calculations in which the absolute value of the energy, being 
a negative power, decreases monotonically to as the relevant geometrical pa- 
rameter increases. In (I77j) the PV force is positive although the PV energy is 
negative; the individual terms in the PD force can have either sign, although 
their energies are all negative; the HD force is zero because the positive HD 
energy is independent of a. Piston analyses center on the cancellation of the 
positive PV force by an external force, since the sum of the other three terms 
can be shown to be negative (see Sec. HJ). 

As a first step toward less naive renormalization, one can keep the cutoff t finite 
and retain the cutoff-dependent terms in (1T6]) . Then the VP term produces a 
force 

f,vP = +5 ^, (78) 

and the terms in (ITT]) are modified in ways that can become significant when 
a or b is not large compared to t. 

Another way to calculate the force is to integrate (Tu(a,y)) over the side of 
the box. We take a moment to verify that the methods are consistent, using 
the appropriate formulas from Sec. 12.31 



Periodic paths: Multiplying (1481) by b to perform the trivial integration, we 
get (after t — > 0) 

b -2{2ja) 2 + {2kb) 2 
F ' k ~ 4^ [(2ja) 2 + (2kb) 2 f 2 ' 



Here j and k are not both ; the terms with one, the other, or neither add 
up to the first three terms in ( 1TT1) . as expected. 

Corner paths: From the energy calculation we know that these terms should 
be zero. Also, from (I5T1) we have T n = unless /3 = £ — \ is nonzero. For the 
f3 term (|52|) . note that the integrand has the form of a total derivative, 

K-2(2bk-2y) 2 d y-bk 



[K + {2bk - 2y) 2 fl 2 dy [K + [2bk - 2y) 2 f/ 2 ' 
Setting x = a and integrating over y, one gets 

F = _P{ (l-*0fr 

tQ]k 7i \ [t 2 + 4(j - l) 2 a 2 + 4b 2 (k - l) 2 ] 3 / 2 
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Fig. 5. Two regions where divergent surface energy appears. 

HQ6 

[t 2 + 4(j - l) 2 a 2 + 4b 2 k 2 } 3 / 2 
The sum of ( 1801) over k from -oo to oo telescopes to 0. 

Side paths: The class of paths that bounce off the side in question (along with 
an even number of additional reflections) have T n identically zero fl57|) . This 
is so even though the shortest such paths (those with k — 0, j — 1) give rise 
to a divergent energy in the region marked (3 in Fig. This matches the 
in (E7J). 

Much more interesting are the paths that bounce off the horizontal walls. From 

dMD-§5D, 




(T n ) t Vjk = ~-rlt 2 + (2j«) 2 + (2kb - 2yf]- b/2 [t 2 ~ 2(2ja) 2 + {2kb - 2y) 2 ] 
+ -[t 2 + (2ja) 2 + {2kb - 2y) 2 ]- 5/2 [t 2 + (2ja) 2 - 2(2kb - 2yf\. 

7T 

For fixed j, k the f3 term is just like the corresponding corner term with j — 1 
replaced by j and the sign changed. Therefore, these two classes of (3 terms 
would cancel when summed over j , even if they did not vanish when summed 
over k as we just saw. 

It remains to integrate the other part of (Tu) t vjk over y from to b. The terms 
with j = lead to a clone of the calculation following ( |671) . In particular, those 
terms for which also k = or 1 are divergent when t — > 0. This divergent 
pressure clearly corresponds, in the case k — 1, to the divergent energy in 
region a associated with paths VP that bounce perpendicularly off the top 
boundary. (From k = comes a corresponding effect at the bottom boundary, 
not indicated in Fig. [5j) That energy is proportional to the length of the box 
and hence gives a force ( 17H1) upon differentiation. 

Finally, one wants to integrate the terms with j ^ and see that they re- 
produce the remaining (VD) term in (1771) . The integral of each term is, at 
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b[32a 4 j 4 +16a 2 b 2 j 2 (k-l) 2 ](k-l) _ bfc[32a 4 j 4 +16a 2 b 2 j 2 fc 2 ] 
J_ [4a 2 j 2 +4b 2 (fc-l) 2 ] 3 / 2 [4a 2 j 2 +4b 2 fc 2 ] 3 / 2 

4^ 16a 4 f 



At first glance it may seem that this expression sums over k to zero, by the 
same telescoping argument used elsewhere. However, unlike those previous 
sums, in this case the individual terms do not approach as \k\ — > oo; rather, 
they go to l/(327ra 2 j 2 ). Taking account of both signs of k and j, one gets the 
force to be 



1 VD 



4E 



C(2) vr 



{ 327ra 2 j 2 8vra 2 48a 2 



as needed. 

Although this exercise may appear redundant, it has underscored two impor- 
tant points. First, doing the calculation in terms of pressure instead of energy 
by no means eliminates the problem of divergences. Second, the divergent 
pressure on a given wall is not associated with the divergent energy adjacent 
to the wall (in region [3 in Fig. EJ. Rather, it goes with the divergent energy 
adjacent to the intersecting perpendicular walls (such as in region a). 



4 The Casimir piston 



The physical significance of the forces calculated in [TJ, in our Sec. I3.4I and in 
much intervening literature has been called into question. For one thing, unlike 
the celebrated sphere calculations of Boyer [52] and others, these calculations 
are unable to take into account any forces coming from the region outside the 
box. In addition, within the framework of ultraviolet-cutoff regularization the 
uncompensated divergent energy proportional to the surface area cannot be 
easily dismissed in deducing the force conjugate to a dimension whose variation 
changes the surface area. In our case the offending energy is that localized in 
the region a in Fig. [5], which is proportional to the length of the box, and the 
corresponding pressure was also observed in Sec. I3.4I in the direct calculation 
of (Tn) on the movable side of the box. 

Cavalcanti [27] proposed to avoid both problems by considering a different sit- 
uation, the piston (Fig. EJ). The interior partition is free to move horizontally, 
and one is to calculate the force upon it. L is to be taken very large com- 
pared to a and 6. The argument now is that the exterior of the apparatus is 
unchanging and hence irrelevant to the force, whereas both interior chambers 
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« L — a 

Fig. 6. A rectangular piston in dimension 2. Its "shaft" has length L — a, effectively 
infinite. The word "piston" refers both to the movable plate at x = a and to the 
model ELS cl whole. 

can be treated exactly. Furthermore, the total of the interior side lengths is 
independent of the piston position, a, so that the surface divergences cancel 
in the calculation of the force. 



Generalizations and variations of this model have been extensively studied 
[2^112^11301153115^^ 

The piston model is not without its own physical problems, because interac- 
tions between the piston plate and the horizontal sides have been ignored. In 
a realistic experiment, because of ordinary Casimir attraction the plate would 
be unstable to striking the tube wall on edge, after which it would collapse 
against one of the walls of the tube. It may be argued that this objection is 
irrelevant to the question of principle that the piston model is designed to 
address; the only degree of freedom one is varying is a, so it is legitimate 
to imagine that the plate is constrained from moving in any other degree of 
freedom. There still exists a Casimir force between the plate and the nearest 
wall, though it is somehow prevented from causing motion. However, one can 
argue by symmetry that this force has no significant horizontal component, 
so that the piston theorists are justified in ignoring it. Nevertheless, in a real 
apparatus there would surely be some friction with the walls, so the feasibility 
of an experiment to verify the piston analysis is questionable. 

Putting these doubts aside, we summarize and recast the Cavalcanti analysis 
in our framework of closed paths. The finite part of the force on the piston 
from the chamber on the left has been calculated in fl77j) . The force coming 
from the shaft on the right can be found from the same formula, with the 
sign reversed, a replaced by L — a, and L taken to infinity; the only term that 
survives is the PV term, 



It exactly cancels the corresponding term in (1771) . leaving PH, PD, and VD 
terms: 

C(3)6 b_ ~ k 2 b 2 -2j 2 a 2 7T 
pist ~ ~ 8vra 3 + 8^ {fa 2 + k 2 b 2 f/ 2 + 48^ ' ( > 
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Here there is no "naive renormalization" as in (1771) . since the divergences (in 
particular, the VP terms) would explicitly cancel if the calculation were done 
for the complete piston before removing the cutoff. 

Cavalcanti [27] rendered fl82|) more illuminating by subjecting it to further 
analysis. If one refrains from the ((3) simplification, the complete sum over 
periodic paths in the (t = 0) energy, flBol) . is 

(% = 4 £ (jV + W)- 3/2 - (83) 

(j,fc)/(0,0) 

From this one can derive two complementary formulas, useful in the respective 
regimes a>6 and a ^ b. (Unfortunately, none of the three formulas for F P i St 
is completely transparent for a ~ b.) 

In the first case, for j — one evaluates the k sum to the term PV (the first 
term in fl66l) ). as before, but for fixed j ^ 0, one applies a known relation 
between the k sum (which is an Epstein zeta function) and a series of mod- 
ified Bessel functions. (This theorem traces back ultimately to the Poisson 
summation formula; see the appendices of [2] and [8].) Thus the PH and PD 
terms together are replaced by the energy terms 




(The individual terms in the sum cannot be associated with individual periodic 
orbits, nor with individual eigenvalues.) Remarkably, the first term of this 
expression precisely cancels the VD term, so that PH, PD, and VD all together 
reduce to the energy 




Since HD does not contribute to the force and the PV force is still cancelled 
by the force from the shaft, the force on the piston is [27, (11)] 




It follows that the piston force (a) is always negative, (b) vanishes exponen- 
tially fast for a b, in contrast to the usual power-law decay of the Casimir 
force. 

Alternatively, one can apply the Epstein-to-Bessel transformation to the j sum 
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for fixed k. That is, PV and PD get replaced by 



486 



7T 




) 



The first term cancels HD (which doesn't contribute to the force anyway); the 
VD term remains (as does PH); and the PV term has been absorbed, so that 
the force from outside the piston is now uncompensated. Thus the total force 
on the piston is J27J (14)] 



As Cavalcanti explains, this form is nicely adapted to understanding the 
regime a < 6, where the standard Casimir result is, of course, recovered in 
the limit. 

Energy is fungible, so one must beware of attributing too fundamental a con- 
nection between particular classes of paths and the observable net forces. The 
striking thing is that the calculations reveal several exact cancellations, not 
all of which can be implemented at the same time. 



5 The Casimir pistol 

The Casimir piston has proved to be a highly illuminating model, but it does 
not settle the issue of the true physical significance of the purely internal 
vacuum pressure on the side of a rectangular cavity. In the piston both the 
divergent (VP) and the positive finite (PV) internal pressure are exactly bal- 
anced by the precisely analogous pressures in the long shaft on the other side 
of the movable plate. This observation does not tell us what would happen if 
the external shaft were not there. The problem of real interest is a rectangular 
box with one side free to move, as indicated schematically in Fig. 0and more 
realistically in Fig. [7J The main question is whether the force on the movable 
side is attractive or repulsive. This is a question about disjoint rigid bodies, so 
it is a meaningful physical problem, just like the piston. (The 3D electromag- 
netic case should be qualitatively similar to the 2D scalar problem.) Another 
urgent question is what happens to the VP divergence now that there are no 
VP paths in the shaft to compensate it; the previous paradox (Sec. [3]) of an 
apparent infinite pressure has reappeared. 

The problem is difficult because there is no reliable analytical calculation of 
the forces acting from outside the box and inside the tiny gaps at the ends of 
the lid. If we momentarily ignore the gaps, it seems unlikely that the external 
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Fig. 7. A rectangular box with one side ("lid") free to move. The box has walls of 
finite width and a finite gap between the lid and the box sides. 

forces would be very large (although we find unconvincing Lukosz's attempt 
PQ to prove this fact by appealing to Weyl's theorem). If one thinks in terms of 
closed paths, paths striking the walls perpendicularly will yield only the usual 
surface divergence, so the only possible source of nontrivial external forces 
is the diffractive paths striking the corners. If this diffractive effect is small, 
therefore, one might expect the force to be repulsive when the plate is exactly 
at the mouth of the box. However, if the plate is located significantly inside or 
outside the box, intuition says the opposite: "inside" we are getting into the 
piston regime, whereas "outside" the Casimir attraction between the nearest 
neighboring regions of the two bodies should be dominant. 

A convincing resolution of this apparent paradox presumably requires a serious 
study of the gap region in a less idealized geometry, as in Fig. [7J It is clear that 
what happens around the gap is very complicated, especially when the plate 
is part-in and part-out as in that figure. One should note that the symmetry 
argument used in Sec.Hlto dismiss the forces in the gap is no longer applicable. 

The uncertainty about the external and gap forces is somewhat alleviated if 
we replace the thin lid by a large rectangular object (Fig. [8}. The piston plate 
has now become more like a bullet or artillery shell. The question now is the 
sign of the force for various values of the five dimensions indicated in Fig. [HJ 
Does a Casimir pistol exist? The advantage of this new problem is that the 
corners of the two bodies are not near each other, so there are no short classical 
paths outside the apparatus (as long as neither d nor e — d is small), even if 
diffractive paths are admitted as classical. (One could eliminate diffractive 
paths (in the sense we are using the term) by replacing the barrel and bullet 
by similarly shaped objects with smooth boundaries.) Like all piston authors, 
we continue to consider only horizontal motion (variation of a) and therefore 
ignore the vertical force between the bullet and the shaft of the barrel. 

We now consider the implications of taking the small gap of width c seriously. 
(We speak only of the gap at the top, but obviously the same remarks apply 
to the one at the bottom.) The first (and motivating) observation is that the 
total side length of the system is now fixed, and hence so is the transverse 
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Fig. 8. The Casimir pistol, consisting of two disjoint, perfectly conducting bodies, 
the barrel and the bullet. It is shown schematically, as in Fig. [51 but the barrel can 
be thought of as having finite thickness, as in Fig. 

extent of the infinite (or cutoff-dependent) surface energy. In particular, the 
energy associated with what we call VP paths (including those striking the 
exterior of the apparatus) is independent of a. The associated paradox is 
thereby removed! 

Our joy in this victory should be short-lived, however. If we take VP paths 
inside the gap seriously, then for consistency we must also take PV paths 
across the gap seriously, and we shall see that they present a serious problem. 

Nonperpendicular paths inside the gap eventually escape from it, so they are 
not "short" and probably can be neglected. (In the model of [33] , all such paths 
escape to infinity and hence can never be closed. In our case it is possible, but 
rare, for such a path to bounce off the left side of the rectangle and return to 
its starting point.) From the point of view of a point r inside the rectangular 
region of area ab, the box now has small "leaks" of width c, but one would not 
expect that to affect its internal Casimir energy significantly. This observation 
could be made quantitative by imitating a calculation in [33], but we shall not 
do so here, because we are interested only in the limit of very small c. 

On the other hand, because the gap c is much smaller than the box dimen- 
sions, a and b, the Casimir energy associated with the rectangle of area cd is 
much greater than that of the box. The principal force associated with this 
gap rectangle is the vertical Casimir attraction between the bullet and the 
barrel, but we have agreed to impose a constraint that makes it irrelevant. 
However, it is the very essence of the piston argument, especially as developed 
by Hertzberg et al. [29 30J, that the proportionality of the Casimir energy to d 
produces a horizontal force, independent of d but proportional to 1/c 2 . (This 
energy is precisely the contribution of the vertical periodic (PV) paths.) In the 
present scenario this force has sign opposite to the Lukosz force in the box, 
because d increases when a decreases, and a larger magnitude than the anal- 
ogous force in the piston scenario, because c < b. Therefore, if we accept all 
the approximations involved in this argument, we are forced to the conclusion 
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that the bullet is sucked into the barrel, not expelled from it. 
Let us list those assumptions. 

(1) The gap between the bullet and the barrel does not significantly affect 
the Lukosz force from the empty part of the pistol chamber. 

(2) There are no significant forces from outside the pistol. 

(3) The effect of the gap can be estimated by ignoring nonperpendicular 
paths and treating the perpendicular paths as usual, as if we had simply 
a pair of parallel plates there. 

Obviously, a trustworthy treatment of this system requires either a numerical 
analysis (for example, by the method of Gies et al. [63]) or, better still, a 
breakthough in the analytical treatment of convex corners. As a positive result 
in this direction, we report that it is possible to compute exact forces and 
torques between bodies of arbitrary shape in weak coupling (for example, 
materials with dielectric constant nearly unity). For example, two thin parallel 
plates of finite length experience an attractive lateral force that tends to cause 
the plates to move to a configuration where they are centered on each other. 
This is the attractive force that tends to increase the length of the gap in 
the pistol. (A recent independent investigation [HI] likewise shows a system 
maximizing the length of a small gap between flat surfaces.) Moreover, in 
addition to the attractive force between the plates, there is a torque exerted 
on one thin plate above a larger plate which tends to cause a rotation of the 
smaller plate about its center of mass so as to favor perpendicular orientation. 
Details have been reported in [39,40] . Since these qualitative conclusions are 
essentially geometrical, they should also hold for strong coupling (Dirichlet 
boundary conditions). 

Finally, let us try to confirm the foregoing conclusions by looking at pressure 
integrals. In principle, one can find the total force on the bullet by integrating 
the appropriate components of (T^J) over the surface of the bullet (or even 
some larger surrounding surface [3TH32] ) . The integral over the back side of the 
bullet is essentially the same as in Sec. [3] (apart from the "infinite" term). On 
the top and bottom sides, the relevant component is T 12 , and a check of the 
formulas (14"9"1) . (15"Ti) . ([SSD? (|BTI) shows that the contributions all vanish. So, one 
would conclude that the pistol fires after all! We believe that the resolution 
of this new paradox is that the crude approximations listed above, although 
they may be permissible for the energy calculation, are simply wrong for the 
pressure calculation. In particular, if the Casimir (or the Lukosz) formulas were 
accurate over the entire gap rectangle, there would be finite jumps at the end 
surfaces of the gap in (Tn) (which is constant and large in the gap, constant 
and smaller in the chamber, and zero in the exterior, in our approximations). 
By the conservation law, there is then a delta function (of x — (endpoint)) 
in (dT 12 /dy). A more realistic calculation would smear out this singularity, 
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probably creating a lump of (T12) that decreases more or less linearly in y 
away from the back corner of the bullet and also downward away from the 
front edge of the barrel. These stress terms would create horizontal forces. 
They are very much like the stresses found in [31, Figures 4(d,e,f)]. 



6 The Casimir pistol with cutoff 

6. 1 Parallel plates revisited 

Our discussion so far has concerned the 2D scalar analogue of the idealized 
perfect-conductor model of the interaction of the electromagnetic field with 
metal bodies. It is generally agreed that the divergences (except for the uni- 
versal volume divergence) encountered in such calculations are the fault of the 
physical failure of that model at high frequencies — equivalently, at length 
scales so small that the material cannot be modeled as a continuum. It is also 
now agreed that the energy divergences, or the corresponding cutoff-dependent 
terms in a calculation with a cutoff, being independent of the bodies' positions, 
do not appear in the forces between rigid conducting bodies. It is sometimes 
forgotten that the idealized Casimir theory runs into physical trouble already 
for rigid bodies, even the canonical scenario of parallel flat plates, when the 
distances become too small. It predicts an energy per unit cross section, £ , pro- 
portional to — aT d for plates with separation a in <i-dimensional space. If taken 
literally, this says, implausibly, that S becomes (negatively) infinite when a 
goes to zero. One would expect instead that in that limit 8 approaches a con- 
stant, since then the space between the plates has disappeared and space is 
filled by the perfectly conducting material. (In fact, the constant should make 
the total energy turn out to be when suitably defined surface energies are 
also taken into account.) 

Barton [HS] has done extensive calculations for dielectric bodies with a po- 
larizability small enough to be treated perturbatively (the opposite regime 
from perfect conductivity). He showed (see also [66]) that a spatial cutoff at 
atomic distances serves to cure the divergences (which otherwise remain even 
in the usual model of quadratic falloff of dielectric constant with frequency 
- e.g., [6"T1I6"%] ). Roughly speaking, the mathematical effect of such a cutoff is 
similar to that of a very rapid, such as exponential, cutoff at high frequency. 
In Barton's theory the total energy per unit cross section does approach as 
a — > when the surface energy is included. In this model, as a — > there is 
a constant attractive force proportional to the energy density of the uniform 
medium, no matter how the latter is regulated |69j . 

More recently, Barton [TUfTT] has developed a plasma model that is more 
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pertinent to the limit of perfect conductivity. It also involves an atomic-scale 
cutoff, but one affecting only the wavelengths parallel to the boundaries. 

Our aim here is to stay in the highly conducting regime and to see whether 
keeping the exponential cutoff parameter t^ 1 finite, at some value typical of 
atomic separations, yields a physically plausible (and divergence-free) model of 
Casimir phenomena. Although ultimately no substitute for serious microscopic 
modeling of conductive materials (an unavoidably nonlinear problem), this 
approach offers hope of rescuing the huge investment that has been made 
into treating vacuum problems (relatively easily) by spectral analysis of linear 
partial differential operators. It also provides a route to understanding the 
gravitational significance of "divergent" local energies and stresses |42j . This 
cutoff should be regarded as analogous to the ad hoc repulsive core in the 
Lennard- Jones potential in atomic physics. A more accurate potential should 
be based on the electronic structure of the atoms; but one would not then 
apply such a potential to, say, nucleon-nucleon scattering. Similarly, a detailed 
theory of real metals is not relevant to hadron bags, cosmological branes, 
thermal fluctuations in soft-matter physics, and other systems where Casimir- 
like effects have been studied. Within our two-dimensional scalar model (which 
is pertinent to all these contexts, if to any) the simple exponential cutoff has 
the advantage of being universal, but we and readers must remain conscious 
that its relevance at small distances to any particular real physical system is 
qualitative at best. We stress again that this atomic-scale cutoff must not be 
confused with the well known decrease of dielectric constant with frequency 
above the plasma frequency; we make no attempt to model the latter, which 
is specific to the electromagnetic scenario. 

In the context of the Casimir pistol, the idea is that the small gap c surrounding 
the bullet must be in the sub-Casimir regime if the other dimensions (a, b, 
d) are in the regime where Casimir effects are significant, and, therefore, the 
deduction in Sec. Oof a dominant attractive force originating in the gap goes 
outside the regime of validity of the theory. Although the cutoff theory has 
no fundamental physical justification, it is probably a bit closer to the truth. 
(Unfortunately, we shall see that no robust conclusion is attainable by this 
route.) 

Maclay and Villarreal [72] proposed this same kind of cutoff and hence ob- 
tained formulas and graphs rather similar to ours in this section. However, they 
identified t with the reciprocal of the plasma frequency rather than, as we do, 
the interatomic spacing, which is typically 100 times smaller (again cf. [65J). 
Other authors [73171] considered (for refutation) an even bigger exponential 
cutoff length, adequate to make the unrenormalized vacuum energy of empty 
space consistent with cosmological observations, and showed that such theo- 
ries predict Casimir repulsion at distances large enough to be refuted by the 
existing laboratory experiments. 
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Fig. 9. Graph of F(r) as a function of r. 

Before studying the pistol, let us look at the attraction between parallel plates 
with the cutoff retained. This could be done easily in any dimension, but for 
coherence in this paper we retain dimension 2. If the separation between plates 
is a, we take the 5 > a limit of (l65j) . in which only the perpendicular paths 
(k = 0) contribute, and divide by b to get energy per unit length: 

c _q- t 2 -2j 2 a 2 

7r^(t 2 + 4j 2 a 2 ) 5 / 2 ' 1 J 



This function behaves in keeping with the idea that the energy or the force 
should be damped when a is comparable to the nanoscale (interatomic spacing) 
represented by t. Of course, when a>t, the effect of t is negligible and (I86I) 
gives the standard result. It is convenient to measure a in units of t. If a — st, 
then irt 2 £ = F(s), where 

00 1 _ 9?' 2 r 2 

F(r)5r £ d + W» (87) 



(see Fig. [9]). F(r) has a zero at ro ~ 0.5888. It has a minimum (a zero of the 
force) at n « 1.0105, with F(n) « -0.02821. At large r, F(r) ~ -C(3)/16r 2 
as in the theory without cutoff. For small r the Euler-Maclaurin formula [751 
(23.1.30)] shows that F{r) ~ | - | + 0(r w ) for arbitrarily large JV. Thus F(0) 
precisely cancels the surface energy from (ITUj) (where P = 2 because we are 
looking at unit cross section on two plates), so that the total energy at s = 
is indeed 0. (But this result may be an accident. It does not happen for the 
Neumann boundary condition. Also, as we shall now observe, s = seems to 
represent material under compression, not a solid block of ordinary material.) 

At r < r\ this model predicts a repulsion. Therefore, it must violate the 
hypotheses of the theorems stating that vacuum forces between symmetrical 
bodies separated by a plane are always attractive [TSfTT] . The argument of 
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Kenneth and Klich [76] refers to the standard dielectric model of the media, 
or a scalar analog thereof, into which our cutoff does not fit. The mathematical 
reason why the theorem of Bachas [77] doesn't apply is less clear, but the key 
physical point is clear from that author's remarks (p. 9094) that a "quantized 
particle does not, strictly speaking, live in one side of the reflecting plane," and 
that the theorem would apply at the quantum level only if (in the terms of our 
scenario) one of the slabs were made of antimatter. The repulsion occurs only 
at separations of the order of the interatomic spacing. Thus the model mocks 
up a more realistic theory in which the two slabs are not cleanly separated, in 
accordance with Barton's remark [551 p. 4088], "[A sharp short-distance cut- 
off,] though a fiction, is a convenient shortcut to somewhere near the truth. At 
small separations, overlap between the electron clouds makes the interatomic 
potential highly repulsive...." We recall also that Ford and Svaiter [21] found 
that a similar effect was induced by a stochastic uncertainty in the position of 
the conducting boundaries, which must in general lead to some probability of 
interpenetration. In short, nobody should be surprised to encounter a repul- 
sion when pushing two slabs of material together. A normal, stable material 
must resist compression. Of course, such repulsion is not a "Casimir effect" ; a 
quantitative study would require detailed modeling of the material, and it is 
a surprise and probably an accident that our crude field-theoretic model gives 
such plausible results in this regime. In particular, the fact that our potential 
minimum occurs at a small positive separation, rather than zero or negative, 
is not to be taken too seriously. 

6.2 Energy in the pistol 

Now we do the energy accounting for the pistol, under the three assumptions 
listed in Sec. [51 (The notation is a slight simplification of that in Sec. [3j) 

Energy in the chamber: According to assumption (1), the contribution Ej> of 
periodic paths is still given by (I65I) with rj = 0. Corner paths can be ignored 
because they make no contribution to the total energy. HP and HD paths 
can be ignored here because they make no contribution to the relevant force 
(their energies being independent of a) . The contribution of VP paths is the 
a-dependent term of flTOj) : 



The contribution of VD paths is given by the generalization of (J7TT) to finite 
cutoff, 
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Energy in the barrel: In the formulas above, (a, b) must be replaced by (d, c), 
and we must multiply by 2 to count both top and bottom gaps. In accordance 
with assumption (3), only vertical paths (j = 0) will be considered. The PV 
paths give 

_2cd~ t 2 -2fc 2 c 2 



The VP paths give 

£VP. = -^, (91) 



of which half belongs to the barrel and half to the bullet. As expected, the 
barrel part combines with 



-Evp + TT-^VP 



1 a + d 



2 v " 4vrf 2 



to yield something independent of a, because a + d is constant. Similarly, the 
bullet part of f[9~Tl) combines with the surface energy of the part of the bullet 
outside the barrel. 

Summary of pistol energy: The only energy terms that contribute to the force 
(under our assumptions) are 

E = Ep + -Evd "i" Ep> (92) 



as listed above. We could differentiate with respect to —a (using = — 
to get the force. All the sums encountered can be expressed in terms of inho- 
mogeneous Epstein zeta functions [6.8J. However, for our purposes it is better 
to analyze the various terms qualitatively. (Quantitatively, we claim nothing 
for the model at short distances anyway.) 



6.3 Asymptotics and numerics for the pistol 

Let c = rt, a = st, b = ut, d = L — a = (/ — s)t (see Fig. [TUl) . We want to 
examine E as a function of s, with r of order unity and s, u, l — s much larger. 
From (!92l) and (165]) we have 

E = Epy + EpH + Eppj + Eypi + Epi 
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Fig. 10. Pistol dimensions in units of t. 

__us™ 1 - 2k 2 u 2 us™ 1-2jV 
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2ms l-2j 2 s 2 -2£; 2 w 2 



A — 1 7 1 \ - 1 - 



vrt £J £i (1 + 4j 2 s 2 + 4Pm 2 ) 5 /2 

g ~ -l+4j 2 s 2 2r(/-s) ~ l-2fcV 
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Let Ep// denote the part of Ep> proportional to s. The other term in Epi 
(proportional to I) is independent of s and hence shall be ignored in further 
discussion of the force on the bullet (including Figs. riTTfT2|) . 

The terms Epy and Ep» are linear functions of s, while the other three terms 
are nonlinear. The linear terms dominate the force at large s, and the main 
point of interest is the confrontation of Epn (Casimir energy in the gap) with 
Epy (identified in Sec. 13.41 as the source of the Lukosz repulsive force; it is 
the term that would give an attractive Casimir force between the upper and 
lower walls of the chamber if those were allowed to move). We shall see that 
genericaliy the P" term is dominant. At small s the nonlinear terms dominate 
and collectively give a function qualitatively similar to that in Fig. [91 Two 
cases are exhibited in Figs. [11] and [12j 

In more detail, by approximating the sums by integrals one can show that 
(when the dimensions other than r are ^> 1) 

C(3)« • r C'U 

ntEpn ~ — =>■ attractive force ~ 5- ; (94) 

165^ s d 



C(2) C 

TctEyp, ~ H =^> repulsive force ~ + — ; (95) 

8s s 2 



ntEpD ~ h ( — — for some function h, such that 



uj u 



C(2) C 

s ^> u ^> 1 TxtEpp, ~ =>- attractive force ~ (96) 

8s s 2 
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Fig. 11. Graphs of linear (dashed) and nonlinear (solid) parts of irtE(s) for r = 1, 
u = 100, I = 500. (a) Small s; linear terms are negligible, (b) Large s; linear terms 
dominate and create an attractive force. 

so that Ey D and E PD cancel to leading order in 1/s, and 

u > 8 > 1 TrtEpD ~ - £p + O ( 4") (97) 



so that the PD force vanishes to leading order in 1/u. In the regime u>s the 
PH force dominates the other nonlinear terms. (In the absence of fine tuning, 
it is still smaller than the linear term unless u ^> s 3 .) In the small-t limit 
(i.e., s ^> 1) the PH force reduces to the first term in (|85p . which is simply 
the standard Casimir force between the left chamber wall and the bullet, 
whereas for general s its energy is E Pii = uF(s)/irt — exactly proportional 
to the parallel-plate function in Fig. [9l Finally, for large s (s ^> w 1 / 3 ) all 
the nonlinear forces are small compared to the linear ones, unless the latter 
happen to cancel. 
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Fig. 12. Graphs of linear (dashed) and nonlinear (solid) parts of irtE(s) for r = 0.5, 
u = 100, I = 500. (a) Small s; linear terms are negligible, (b) Large s; linear terms 
dominate and create a repulsive force. 

Indeed, the forces arising from the linear terms are 

nt 2 Fp Y = -F(u), irt 2 Fp„ = +2F(r), (98) 



where F is defined by (1871) . Recall that F has a zero at r ~ 0.6 and a minimum 
at ri ~ 1 and rapidly approaches at large r. When s ^> w 1 / 3 , Fpn exceeds all 
nonlinear forces, in particular the PH force (which, we have seen, dominates 
the nonlinear forces if s u) . Whether the total force is attractive or repulsive 
at large s is determined by the relative size of the two constant forces in (198|) , 
and hence on the value of r, F{u) being small and negative. (1) For r < r u 
and r not too close to r , Fpn dominates and the total force is attractive. This 
is the regime in which the cutoff model seems most trustworthy physically. 
In particular, it contains the point r\ , which one might regard as the most 
"natural" value, corresponding to two blocks of material in relaxed contact, 
their effective surfaces separated by the typical interatomic spacing. (2) For 
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< r < r and r not too close to r$ , Fp» again dominates and the total force 
is repulsive. In particular, if r = the two gap forces (P" and VP') cancel 
and the force can be attributed to the negative surface energy in the region 
a (Fig. E]) and the similar region next to the part of the bullet outside the 
barrel. (3) If r = Tq , Fpn vanishes and the long-range force is purely Fpy, the 
repulsive Lukosz force. This result is what the piston model was designed to 
achieve — a gedankenexperiment showing that the Lukosz result has, at least 
in principle, some physical reality. Unfortunately, that result is attainable only 
by fine-tuning and, moreover, by pushing r into a regime where the physical 
relevance of the cutoff model is questionable. (4) For a special value close to 
tq , namely 

F(u) C(3) , N 



the long-range linear force vanishes. In this scenario the force is the sum of the 
PH, PD, and VD terms; in the small-cutoff regime of interest, it is approxi- 
mately the piston force (1821) . This force is always attractive, but exponentially 
weak at large s. 

In short, the gap plays a spoiler role somewhat like that of the outer shaft 
in the piston model. In the piston, where shaft and chamber have the same 
width, the shaft force precisely cancels the related repulsive force from the VP 
paths in the chamber, leaving the Casimir-like force ( |82|) . The same occurs for 
the pistol in scenario 4, but in the more plausible scenario, 1, the (attractive, 
a-independent) gap force overwhelms the interior VP force, precisely because 
the gap is narrower than the chamber. In any case, the force arising from out- 
side the chamber depends, not surprisingly, on the geometrical configuration 
outside the chamber, while the force arising inside is fixed by the geometry of 
the chamber. 

For the reason mentioned in connection with parallel plates, the fact that the 
gap force becomes repulsive at all at finite gap size may be an artifact of 
the cutoff model. What one can say is that the even cruder model of perfect 
reflection also displays an artifact, in the form of an attractive force that 
diverges as the gap size approaches zero. In a more realistic model taking into 
account the interactions related to condensed matter physics one would expect 
the gap force to be reduced, if never reversed. 



7 Conclusions 



We have presented a thorough analysis of the vacuum expectation value of 
the stress-energy-momentum tensor in a rectangle. The calculational methods 
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involve an exponential ultraviolet cutoff and a sum over images (or closed 
reflecting paths). Here we have treated a two-dimensional scalar field; the 
extension to three dimensions and electromagnetism is straightforward and 
under way. Formulas are presented for all tensor components, T^r), for ar- 
bitrary combinations of Dirichlet and Neumann boundaries, arbitrary values 
of the curvature coupling £, and arbitrary values of the cutoff parameter, in- 
cluding the limit where the cutoff is removed. Forces (which are independent 
of £) have been consistently calculated both by differentiating energy and by 
integrating pressure. 

Studying the local energy density and stresses (rather than just total energy), 
using a physically motivated ultraviolet cutoff (rather than an "analytic" reg- 
ularization scheme), and studying separately the contributions from various 
classes of specularly reflecting paths all help to interpret the physics, espe- 
cially the roles of boundaries and corners. Within a cutoff framework one has 
a clear and consistent definition of energy densities and forces. When differ- 
ent configurations of rigid bodies are compared and all contributions (from 
inside and outside) are included, one always finds a cancellation of the energy 
divergences and hence an unambiguous force in the limit of no cutoff. The 
decomposition by paths helps one to understand better the cancellations of 
divergent terms and often to understand intuitively the sign of the Casimir 
force. Most strikingly, the force on one side of the rectangle includes important 
repulsive components associated with paths parallel to that side: a divergent 
term from short paths that reflect from the perpendicular sides, and a finite, 
constant term from periodic paths between the two perpendicular sides. In 
piston geometries these forces are cancelled by counterpart terms from the 
exterior of the rectangle, but in more general circumstances the problem of 
their physical interpretation must be taken seriously. 

In the later sections of the paper we discuss geometries in which the vacuum 
forces from inside a rectangle might be rigorously exhibited. The box with 
a loose lid (Fig. [7j) is closest to what one wants to understand, but accu- 
rate calculation of the external edge and, especially, corner effects remains 
impractical for now (at least, beyond the scope of the present paper). The 
piston model (Fig. [5]) studied by previous authors is rigorous and exact, but 
it obscures the point at issue by adding an external shaft. Our attempt to 
compromise these two scenarios is the pistol (Fig. [SJ, which unfortunately did 
not yield a robust result. The force on the pistol depends sensitively on the 
cutoff length, as compared to the width of the gap between the bullet and 
the barrel. The only regime in which our quantitative analysis (extrapolated 
to 3D electromagnetism) can be regarded as physically trustworthy is that 
where the gap is small but still larger than the cutoff; there the behavior is 
cutoff-independent but the force is attractive. Scenarios where the net force 
is repulsive (in particular, one where the gap force vanishes) do exist, but re- 
quire entering the regime where the calculations cannot be taken seriously on 
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a quantitative level because one does not know what the correct ultraviolet 
cutoff behavior is (and because stiction and friction are likely to be the dom- 
inant effects there); furthermore, making the gap force zero or small requires 
fine tuning within this regime. Nevertheless, although no quantitative claims 
can be made for our model (pistol + cutoff) in that regime, we do submit that 
the model is closer to the physical truth than either a model without cutoff 
(which would predict infinite energies) or an analytic regularization that hides 
the divergences from the beginning. Furthermore, while the repulsive Lukosz 
component of the force is robust, the force opposing it is dependent on the 
scenario considered (e.g., piston vs. pistol, or wide gap vs. narrow) and could 
in principle be controlled to demonstrate the reality of the Lukosz force, even 
if the net force is attractive in all practical experiments. 
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